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Abstract 


This  article  reports  on  a  positron  computerized  tomography  camera  system 
using  multiv'ire  proportional  chambers  (MWPC)  as  detectors.  This  system  is 
composed  of  two  high-density  MWPC  gamma-ray  detectors,  an  electronic  readout 
system  and  a  computer  foi  data  processing.  Three-dimensional  tomography  is 
obtained.  The  imaging  matrix  is  64  x  64  x  16. 


1.  Introduction 

The  appearance  of  x-ray  computerized  tomography  (XCT)  is  a  great 
pioneering  work  in  the  field  of  medical  imaging,  encouraging  those  involved  in 
nuclear  medicine  to  bring  to  fulfillment  their  aspirations  for  three- 
dimensional  radiography.  Nuclides  emitting  positrons  are  used  as  the  tracer; 
electrons  from  within  the  human  body  and  these  positrons  are  annihilated,  and 
a  pair  of  gamma  rays  is  emitted  in  opposite  directions.  Use  of  the 
appropriate  technique  to  measure  these  two  rays  makes  it  possible  to  obtain 
information  on  direction  and  position.  Tomographs  produced  by  positrons 
(called  PECT)  are  one  of  the  best  means  available  to  the  medical  imaging  field 
for  realizing  three-dimensional  description;  what  they  provide  is  not  merely 
anatomical  information;  their  more  important  value  lies  in  the  transmission  of 
information  on  normal  and  abnormal  human  chemistry  and  physiology  [1].  In 
nuclear  medicine,  the  use  of  PECT  to  investigate  the  circulatory  map  of  the 
brain  and  the  foot,  and  to  investigate  the  metabolic  process  of  the  cerebrum, 
the  heart,  the  liver,  the  thyroid  and  other  organs,  and  especially  to  study 
diseases  of  the  nervous  system  of  the  brain,  has  conspicuous  merits. 


1 


In  the  past  ten  years,  the  development  and  use  of  PECT  has  been  given  the 
greatest  possible  emphasis  by  physicists  and  medical  scientists.  Thousands 
have  devoted  their  efforts  tc  research  and  development  in  the  field.  The 
Berkeley  (U.S.A.)  Laboratory  [2].  the  European  Center  for  Nuclear  Research 
(CERN)  [3],  and  England's  Rutherford  Laboratory  [4]  have  developed 
individually  different  forms  of  raultiwire  positron  computerized  tomography  for 
nuclear  medical  research  [5] .  In  addition,  first  and  second  generation 
multiwire  proportional  chamber  {MWPO  PECT  is  directly  used  in  the  field  of 
solid  body  investigation  in  physics  to  measure  Fermi  surfaces  as  well  as  to 
determine  Che  moirentum,  distribution  and  concentration  of  electrons  in  metals 
and  alloys  [6].  Ve  began  tc  develop  MWPC  PECT  in  1983  and  after  three  years' 
effort  have  achieved  highly  satisfactory  results. 

This  article  describes  in  detail  a  MWPC  positron  camera  system.  It 
includes  a  high-density  multiwire  chamber,  readout  electronics,  computer  on¬ 
line  data  gathering,  and  an  image  reconstruction  process.  Using  a  liquid 
positron  source  to  perform  simulated  static  radiography,  we  achieved  three- 
dimensional  tomography  with  a  clear  image  and  stable  system  operation. 

2.  Detectors 

The  two-dimensional  readout  high-density  MVPC  is  the  core  of  the  positron 
camera.  It  is  composed  of  a  multiwire  chamber  and  a  multilayer  lead-sheet 
gamma-ray  transforming  body  [7] .  The  multiwire  chamber  has  three  wire 
surfaces:  the  central  surface  is  the  anode,  and  the  two  cathode  surfaces  are 

positioned  on  both  sides  of  the  anode.  Reactive  signals  on  the  wires  of  the 
cathode  are  used  tc  record  incident  particles'  X  and  Y  coordinates. 

Because  MVPC  is  a  gas  detector  operating  under  atmospheric  pressure,  it 
requires  an  effective  detection  capacity  equivalent  to  a  511  keV  photon;  it  is 
necessary  first  to  transform  the  photons  into  electrons.  Based  on  this 
consideration,  a  gamma-ray  transforming  body  with  an  effective  surface  area  of 
20  cm  X  20  cm  and  a  thickness  of  6  mm  was  designed.  In  order  to  facilitate 
the  entry  of  photons/electrons  produced  in  different  lead  layers  in  the 
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transforming  body  into  the  sensitive  area  of  the  multiwire  chamber,  it  is 
necessary  to  add  a  drift  electric  field  to  the  transforming  body. 

An  automatic  air-roixing  system  supplies  the  working  gas.  The  gas  is  able 
to  move  evenly  and  without  obstruction  through  the  transforming  body  and 
sensitive  area  of  the  MWPC.  In  order  to  prevent  interference  from  the 
external  environment,  the  detector  has  an  effective  electronic  screen. 


There  are  many  methods  of  reading  out  the  cathode  signal  of  the  MWPC;  in 
order  to  cut  down  the  number  of  amplifiers  and  reduce  the  demands  on  the 
quality  of  the  front-rank  amplifier,  we  have  chosen  the  cathode  wire  coding 
readout  method  [S]  .  Each  cathode  surface  only  needs  a  28-route  amplifier  to 
realize  the  goal  of  reading  out  192  signal  paths. 


The  major  technical  specifications 

Sensitive  area: 

Working  gas: 

Anode  high  voltage: 

Gamma-ray  transforming  body 
drift  voltage: 

Readout  style: 

Spatial  resolution: 

Detection  efficiency  for 
511  keV  photon: 


for  the  detector  are  as  follows; 
20  X  20  cm^ 

80%  Ar  +  20%  CO2 
2.5  kV 

-1.0  kV 

Reactive  pulse  cluster  center 
1  mm 

5.6% 


3.  Electronic  Readout  System 

The  schematic  MWPC  positron  camera  system  is  shown  in  Fig.  1.  The  MWPC's 
100  anode  wires  are  joined  together  and  linked  with  a  voltage  amplifier  [9]  . 
When  the  anode  signal  has  been  amplified,  it  passes  through  a  constant 
comparison  time-setting  filter  device  (CFD)  and  reaches  a  fast  coincidence 
circuit;  the  coincidence  output  signal  passes  out  through  an  exit,  giving  all 
coincidence  registers  a  signal  to  open.  The  192  signal  wires  of  each  cathode 
surface  are  organized  into  four  groups  of  24  pairs,  and  each  group  and  pair  is 
individually  coupled  with  an  electric  charge  sensitive  front-rank  amplifier. 
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Fig.  1.  Schematic  for  MWPC  positron  camera  system. 

Key:  (1)  MWPC  I;  (2)  object;  (3)  MWPC  II;  (4)  marker;  (5)  coincidence;  (6)  xi 
28  routes:  (7)  yi  28  routes;  (8)  delay;  (9)  gate  producer;  (10)  X2  28  routes; 
(11)  y2  28  routes:  (12)  coincidence  register;  (13)  exit;  (14)  coincidence 
register:  (15)  connector;  (16)  LSI-11/23  computer;  (17)  printer;  (18)  terminal. 


Using  an  electric  charge  sensitive  amplifier,  it  is  possible  to  eliminate 
efficiently  the  undesirable  influence  caused  by  each  group/pair's  output 
terminal  stray  capacitance;  at  the  same  time  it  raises  the  sensitivity  for 
picking  up  signals. 

The  front  rank  amplifier's  output  signals  are  sent  through  the  main 
amplifier  and  screening  mechanism  to  a  Lecroy  4448  coincidence  register.  When 
the  gamma  rays  emitted  as  a  pair  by  the  annihilation  of  a  positron  and  an 
electron  are  detected  individually  by  the  detectors,  each  cathode  surface  has 
a  group  of  response  pulse  outputs  which,  with  the  anode  trigger  signal, 
implement  coincidence  registration;  and  further,  passing  through  a  2912  Uni-Q 
communication  line  switching  device  and  a  CAMAC  3912  controller,  read  out  the 
coincidence  register's  data  on  the  internal  storage  of  the  LSI-11/23  computer. 
The  coincidence  register  is  a  single-width  CAMAC  insert  having  a  16  x  3 
register  path.  Recording  the  signals  of  the  MWPC's  four  cathodes  are  a  total 
of  three  4448  registers  occupying  three  positions  in  the  CAMAC  controller. 
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The  computer  uses  the  single-user  RT-11  operating  system,  which  has  a 
rapid  real-time  response  capacity.  The  entire  on-line  transport  program  uses 
the  dependent  program  adjusted  use  method,  the  user  compiling  the  main  program 
in  accordance  with  actual  demand.  The  computer,  on  the  basis  of  the 
et'aluation/selection  criteria  and  the  reactive  pulse  cluster  center  method 
readout  principle,  seeks  out  each  cathode  pulse  cluster's  center,  and 
establishes  the  spatial  coordinates  (xi  ,  yi ) ,  (X2 ,  y2 )  of  an  authentic  photon 
pair  instance. 

4.  Computerized  Image  Reconstruction 

For  image  reconstruction,  we  selected  the  reprojection  wave  filter- 
method.  If  we  take  the  coordinates  (xi  ,  yi  ) ,  (X2 ,  y2 )  of  the  two  gamma  rays 
as  measured  by  the  detector,  join  them  with  a  straight  line,  and  relate  the 
line  to  the  tomographic  plane  of  the  object  under  observation,  the  line  will 
necessarily  pass  through  the  annihilation  point  of  e*  e- ,  called  the  focal 
point.  All  instances  together  comprise  a  set  of  focal  points,  which 
constitute  the  focal  plane,  and  this  is  the  above-mentioned  tomographic 
plane's  reprojected  image.  Obviously,  the  reprojected  image,  in  addition  to 
the  actual  image  formed  by  the  annihilation  points  on  the  plane  under 
consideration,  disc  contains  background  produced  on  the  tomographic  plane  by 
cases  of  annihilation  on  other  planes.  The  process  of  using  mathematics  to 
eliminate  the  background  is  called  the  wave  filter  process;  the  image  after 
undergoing  wave  filtering  reflects  the  distribution  of  atomic  particles  on  the 
tomographic  plane  [10]  .  The  main  steps  in  the  processing  of  the  data  are  as 
follows : 

1)  Using  the  large  quantity  of  measured  data,  establish  the  reprojection 
density  distribution  function  for  each  tomographic  plane;  establish  the 
reprojection  matrix's  Fourier  transformation. 

2)  In  the  transformation  matrix  each  frequency's  amplitude  is  multiplied 
by  the  value  of  the  corresponding  point  response  function,  to  implement  the 
wave  filter  process. 

3)  An  inverse  Fourier  transformation  is  implemented  on  the  filtered 
matrix  to  return  to  actual  space. 
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4)  The  iiTiage  is  smoothed  and  output. 


Under  the  conditions  of  our  experimentation  the  investigated  bodies  were 
placed  between  two  parallel  detectors.  Reprojection  was  a  process  of 
projecting  the  two-dimensional  data  provided  by  the  MWPC's  in  three- 
dimensional  space.  The  mathematical  description  for  realizing  the  reprojected 
image  wave  filter  is  the  Fredholm  formula: 

'  ( 1 ) 
/v(^,  .V,  z)  =  I II  /(*'.>'.*')•  z')dx'dy'dz' 

Here,  I'x',  y,  z)  represents  the  reprojected  image,  that  is  the  experimentally 
measured  radioactivity  distribution:  f(x‘,  y',  z')  represents  the  actual 
three-dimensional  image,  that  is,  the  radioactivity  distribution  that  requires 
solution:  k(x.  y.  z.:  x',  y‘,  z')  is  the  Green  function  for  the  imaging 
system,  also  called  the  point  response  function. 


The  point  response  function  h(x,  y,  z;  x' ,  y',  z')  ordinarily  is  not  a 
spatial  constant.  If  a  sufficient  number  of  limits  arc  added,  for  example 
limitations  to  the  size  of  the  reconstructed  area,  it  causes  the  camera  to  s^t 
the  somewhat  random  reception  angles  in  the  reconstruction  area  approximately 
equal;  large  angle  instances  disappear,  so  that  edge  effect  is  lost,  and  so 
on.  Thus  the  pom:  response  coefficient  becomes  a  spacial  constant,  which  is 
only  the  function  of  /x  -  x').  (y  -  y'),  (z  -  z’),  and  has  no  relation  with 
the  value  of  any  of  the  points  (x,  y,  z:  x' ,  y',  z').  Hereupon  the  integral 
equation  (1)  becomes  the  three-dimensional  bounded  curved  volume  integral: 


/«(•>-,  >,s)  “  I  [  I  y  y  y  z)  ’  —  x  ,  y  —  y  y  z  —  z^dxdy  dz' 

Jo  Jo  Jo 

—  y.  *)  •  A(x,  y,  r). 


[2) 


Using  the  curved  volume  theorem,  formula  (2)  becomes  the  product  of  the 
frequency  range; 

Fu^K.y  Kyy  K .)  «  Kyy  K .)  .  HiK.y  Kyy  K  y) 


Here,  F^.  /V.,  and  H  are  the  Fourier  transformations  ot-^nix,  y,  z),-ffi(x,  y, 
zj,  and  hi'x,  y,  z) ,  respectively. 
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Formula  (3)  can  be  manipulated  to  produce 

FyjK.,  K,,  K.) 


F,(K.,K,,K.) 


HiK.,K,,  K.) 


Using  the  inverse  Fourier  transformation,  we  obtain  the  desired  radioactive 
nuclide  distribution: 


^ »  >1  ' ) 


F t  A.’>,  A’,) 


^  .'J  A,  1/ Ay  J  A . 


(4) 


For  the  two  fixed,  limited  angle  reception  MWPC  position  cameras  used  in 
our  work,  the  point  response  function  is: 


A  ( ^ ,  .V  ,  2  ) 


cosQ 

l%r^ 


X 

ixW  +  y'  +  ' 


(5) 


-  IS  the  angle  of  incidence  of  the  line  of  annihilation  with  z  axis,  cos^  is 
the  ruling  facto:.  The  point  response  function  is  as  shown  in  Fig.  2. 


Fig. 


Point  response  function.  Key: 


(1)  detector  I;  (2)  detector  II. 


Because  the  camera  has  a  limited  angle  of  reception,  the  frequency  space 
point  response  function  K^,  K^)  has  a  frequency  cone  equal  to  zero; 

within  this  cone  the  frequency  cannot  be  measured.  The  reason  lies  in  the 
fact  that  the  camera  loses  the  pair  of  gamma  rays  produced  by  the  annihilation 
of  positrons  outside  the  angle  of  reception.  The  reverse  wave  filter  G(Ki_, 

A'i_,  Ki_)  IS  defined  in  the  broad  sense  as  the  opposite  of  the  wave  filter  H(Ks, 

H-'iK.,  K,,  K.) 

C(,K,  »Ay,A,)“|  ^  ^  Q  (f’) 
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(7) 


Merging  (6)  inti  formula  (3)  we  obtain: 

f  K,,  K,)  “  Fy(.K,,  K,,  K,)  •  G(^K,y  K,,  K,) 

On  the  reverse  wave  filter  G(K^,  Ki_,  K2_) ,  multiplied  by  the  window  function 
V  .Ky_.  A'^.  Kz_  ,  It  IS  possible  effectively  to  control  the  effect  of  the  wave 
filter  magnifying  the  noise  within  the  above-mentioned  frequency  range,  and 
thereby  to  improve  the  quality  of  the  image. 

In  the  process  of  reconstructing  the  image,  the  reverse  wave  filter  with 
a  window  function  is  shown  as: 

G»(A„  Ah,  A.)  -  G(A.,  A,,  A.)  •  H  (A,,  A,,  A.)  (S' 


Sc  forrr.ula  i7)  becomes: 

f  *(  A, ,  A, ,  A. )  “  F A, ,  A,)  •  C^{  K , ,  K , ,  A  , ) 


(9) 


In  this  manner  it  is  possible  to  go  through  an  inverse  Fourier  transfor¬ 
mation,  obtaining  a  post-reconstruction  radioactive  nuclide  distribution: 


/*(r,  v,c) 


A  *( Tv , ,  Ky ,  A,) 


(10) 


The  image  reconstructed  by  the  computer  is  output  clearly  in  shades  of  gray. 

In  order  to  study  results  of  image  reconstruction  and  the  influence  of 
various  parameters  on  image  quality,  we  have  set  up  a  computer  simulated 
mathe.T.atical  pattern  similar  in  shape  to  a  human  cranium.  It  consists  of  the 
Shull,  the  right  and  left  ventricles,  cavities  affected  by  disease,  and 
swollen  lumps.  Each  part  is  provided  with  different  density  distributions  of 
"radioactive  nuclides."  Using  a  computer  we  undertook  a  random  sampling  on 
th<^  m,odel  and  performed  reprojection  and  image  reconstruction  processing.  A 
total  of  10-  "instances"  was  calculated.  Figure  3  shows  the  eighth  and  ninth 
section  images,  including  the  original  image,  the  reprojection,  and  the 
reconstructed  image.  The  results  are  satisfactory. 
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ated  image 


(a)  original;  (b) 


5.  Experimental  Results 


The  MWPC  positron  camera  system  is  shown  in  a  photograph  in  Fig.  4.  It 


Fig.  4.  Multiwire  proportional  chamber  positron  camera  system. 

includes  the  detectors,  high  voltage  electric  supply,  gas  system,  NIM  and 
CAMAC  electronics  with  an  LSI-11/23  computer  having  256  KB  of  memory  and  a  40 
MB  hard  disk,  a  terminal  and  a  printer.  The  two-high  density  MWPC's  are  set 
up  parallel  to  each  other  at  a  distance  of  24  cm.  The  object  to  be  imaged  is 
an  iron-sheet  cylinder  10  cm  in  diameter  and  15  cm  high.  In  the  cylinder  is  a 
columnar  fluid  source  of  positrons,  ^^Na,  with  a  strength  of  200  microcuries, 
enclosed  in  a  glass  bottle  around  which  was  wrapped  a  sheet  of  thin  lead. 

This  kind  of  containment  had  no  great  effect  on  the  penetration  of  the  511  keV 
photon  produced  by  the  e*e-  annihilation. 

The  dimensions  of  the  imaging  space  were  established  as  64  x  64  x  16,  the 
computer  on  line  collected  3  x  10^  instances.  With  data  processing  and  visual 
output,  16  reprojections  and  reconstructed  images  were  produced;  Figure  5 
shows  the  results.  Because  of  the  limitation  on  the  length  of  this  article, 
we  have  selected  only  cuts  5  and  6  of  the  object  center  and  cuts  13  and  14  of 
Its  extreme  end.  For  comparison,  the  reprojections  and  the  reconstructions 
are  shown  in  pairs.  The  three-dimensional  size  of  the  2*Na  positron  source, 
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as  obtained  from  the  reconstructed  image,  is  within  a  3  mm  error  range  in 
comparison  with  the  size  of  the  actual  object. 


/ 


(a) 


(b) 


Fig.  5.  Computerized  reprojection  and  reconstructed  image  for  a  simulated 
object,  (a)  reprojection  for  cuts  5  and  6;  (b)  reconstructed  image  for  cuts 
5  and  6. 
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Fig.  5.  Computerized  reprojections  and  reconstructed  images. 

(c)  reprojection  for  cuts  13  and  14;  (d)  reconstructed  images  for  cuts  13 

and  14. 


Under  the  conditions  of  our  experiment,  the  use  of  the  LSI-11/23  to  make 
a  reprojection  and  reconstructed  image  required  about  one  hour;  the  use  of  a 
VAX-11/780  required  only  two  minutes.  To  increase  the  data  processing  speed 
of  the  LSI-11/23,  it  would  be  necessary  to  equip  it  with  a  high-speed 
calculation  upgrade,  which  could  reduce  the  processing  time  required  to  a 
matter  of  minutes. 

Our  success  in  constructing  images  with  simulated  objects  has  laid  a  good 
foundation  for  trials  with  living  objects.  With  improvements,  the  system  can 
be  used  in  measuring  two  gamma  angle  relations  for  solid  physics 
investigations . 
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